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ABSTRACT 

In  a  new  sonar  bearings-only  solution  method.  Dr.  D.  C.  Bossard 
of  Daniel  Wagner  Associates  achieved  quite  spectacular  reduction  in 
i'ange  errors  on  a  zigging  target,  one-sixth  those  of  the  usual 
(unweighted)  CHURN  method.  His  method  yields  time-corrected  range 
(value  at  time  when  expected  error  is  least)  and  weights  observations 
according  to  assumed  zig  statistics.  Bossard  also  advocates  extrapo¬ 
lating  favorably-chosen  time-corrected  ranges  to  obtain  present 
range. 

We  find  that  the  CHURN,  with  weights  equivalent  to  Bossard 's, 
achieves  equally  small  time-corrected  range  errors,  and  errors  at 
solution  time  one- third  those  of  the  usual  CHURN.  Random  bearing 
noise,  however,  seriously  degrades  solutions  using  Bossard  weights, 
even  without  zigs,  in  which  case  the  unweighted  CHURN  is  optimum. 

For  combinations  of  zigs  and  bearing  noise,  optimum  combined  weighting 
functions  exist. 

Unsuccessful  attempts  were  made  to  use  data  available  to  the 
tracking  ship,  e.g.  autocorrelation  of  solution  residuals,  for 
selecting  optimum  weighting.  Autocorrelation  was  also  probed  for 
zig  detection  clues  without  success. 

Results  obtained  by  extrapolating  pairs  of  time-corrected 
range  to  present  time  were  about  equally  as  good  as  from  single 
solutions  using  the  same  data. 
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We  conclude  that  Bossard's  important  contribution  is  to  show 
the  effectiveness  of  appropriate  statistical  weighting. 

Further  efforts  should  deal  with  on-board  methods  for  optimizing 
weights,  and  the  benefits  of  weighting  for  other  error  sources 
(bearing  bias,  own  ship  position).  The  results  would  be  applicable 
not  only  to  the  CHURN,  but  also  to  the  newer  optimal  filter  methods. 
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I.  INTRODUCTION 

The  proposal  by  Wagner  Associates  to  the  Naval  Ordnance  Systems 
Command  (0RD-0521B)  ,  reference  2,  has*- 1  on  the  work  of  Dr.  D.  C.  Bossard, 
concerns  passive  bearings-only  ranging  on  maneuvering  targets.  An 
example  is  given,  using  synthetic  noise-free  data,  in  vdiich  Bossard 's 
solution  method  shows  range  errors  at  best  range  time  only  one-sixth 
as  large  as  the  corresponding  errors  of  the  CHURN  method  ,s  generally 
implemented . 

We  have  verified  Bossard 's  results  by  duplication.  In  addition 
we  have  obtained  equally  good  results  with  the  CHURN  modified  by  the 
application  of  statistical  weighting  equivalent  to  that  used  by 
Bossard.  The  conclusion  is  therefore  drawn  that  the  power  of  Bossard 's 
method  lies,  not  in  his  novel  analysis  using  bearing  quadruples,  but 
in  weighting  observations  appropriately  for  the  process  (maneuver  in 
this  case)  which  generates  the  errors. 

In  the  Bossard  method  a  moving  window  of  fixed  time  span  is  used, 

A  four -bearing  range  equation  is  derived  which  includes  an  explicit 
error  term.  The  derivation  continues  by  formally  summing  all  such 
equations  obtained  from  suitable  sets  of  bearing  quadruples  drawn 
from  the  window,  weighting  each  addend  in  accordance  with  the 
covariance  of  its  error  term.  Finally,  the  equation  obtained  by 
summing  is  modified  by  adjusting  the  time  parameter  so  that  th.se 
terms  involving  6^1  (target  speed  in  the  line  of  sight)  offset 
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each  other.  Thus  an  equation  for  the  time-corrected  range  for  the 
window  is  obtained. 

The  error  term  appearing  in  Bossard's  four-beari.ng  range 

* 

equation  is  essentially  a  function  of  E,  the  time  derivative  of 
the  residual  E  dealt  with  in  the  CfflJRN  (E  is  cross-range  miss 
distance  of  observed  bearing  lino  from  estimated  target  position) . 

As  discussed  more  fully  later  in  this  report,  the  zig  strategy 
assumed  in  the  geometry  used  here  and  also  by  Bossard  in  reference  2 
generates  an  exponential  covariance  for  E.  Weighting  appropriate 
to  this  covariance  was  used  by  Bossard  to  obtain  his  favorable 
results,  and  by  us  in  duplicating  his  results.  Appendix  A  contains 
the  computing  algorithm  used  at  the  Naval  Weapons  Laboratory  for 
this  purpose. 

We  have  derived  the  covariance  for  E  which  corresponds  to  the 
same  zig  strategy  (see  part  III).  This  covariance  turns  out  to  be 
non-stationary  and  correlated,  therefore  much  different  from  the 
stationary,  uncorrelated  statistics  usually  assumed  in  applying 
CHURN.  Use  of  the  appropriate  covariance  in  CHURN  makes  the  time- 
corrected  range  solutions  agree  rather  closely  with  those  of  Bossard's 
method. 
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Weighting  in  accordance  with  the  covariance  generated  by  a 
zigging  target  evidently  has  some  value  even  for  certain  quite 
different  error  distributions,  as  for  example,  from  the  sinuous 
target  motion  also  tested.  When  the  source  of  errors  includes 
uncorrelated  bearing  noise,  however,  solutions  are  badly  degraded. 
II.  WEIGHTED  GHURW  AND  HME-CORRECTED  RANGE 

The  CHURN  method  estimates  four  target  parameters 
and  a^,  east  and  north  components  of  speed 
a2  and  a^,  east  and  north  coordinates  at  t  =  0 
which  are  assumed  to  be  constant  during  the  solution  data  span. 

From  these  four  parameters,  along  with  own  ship  motion,  all  of  the 
familiar  target  parameters  such  as  range  are  easily  derived.  The 
CHURN  solution  minimizes  the  weighted  sum  of  squares  of  the  residual 

+  a2)cos  -  (a3tj|^  +  a^)sin  B4  -  B^  +  y^j^sin  B^^ 

Of  symbols  appearing  here  and  in  Figure  1, 

Rjj^(a)  is  range  at  time  t^  computed  from  a-^  ...  a^  and  Xq^^,  y^^ 
B^  is  observed  bearing  at  time  t^ 

Xoi,  ygi  are  own  ship  coordinates  at  time  t^. 

Inspection  of  Figure  1  shows  that  Ej;  is  the  miss  distance  of  the 
observed  bearing  line  from  the  computed  target  position.  The  set 
of  equations  for  all  of  the  Ej^  (i  =  1,  2  ...  n)  can  be  represented 
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in  matrix  form 


fl" 

tj^cosB], 

cosB]^  -t^sinBj^  -sinBj^ 

®1 

a<s 

XoiCOsBi-yoiSinBi 

Ei 

= 

tj^cosBi 

cosB.  -t- sinB.  -sinB,- 
111  i 

^4 

- 

XoiCOsBi-yoiSinBi 

En 

tnCOsBn 

cosBn  -tj^sinBj^  -sinB^ 

WosBn-yon^i^Bn 

and,  with  obvious  substitutions,  in  the  more  compact  matrix  notation 

E  =  Aa  -  d 

For  minimum  variance-covariance  of  a,  two  conditions  need  to 
be  met: 

T 

E  WE  =  minimum  over  a 
W  =  constant  x  (G^) 

where  is  the  covariance  of  E.  Equations  expressing  the  first 
condition  are  obtained  by  setting 

-2-  (e’^WE)  =0  1  =  1,  2,  3,  4 

a.i 

then  substituting  for  E,  and  carrying  out  the  differentiation. 

The  re'iulting  four  normal  equations  can  be  written 

=  g 

where  i/'  *  a'^WA  and  g  =  A^Wd 

If  the  residuals  E  are  stationary  and  uncorrelated,  as  has 
usually  been  assumed  in  previous  applications  of  the  CHURN,  then 
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Using  this  relation  to  expand  the  expression  to  be  minimized,  we 
obtain 

t^(  +  2t(  *^44^  minimum  over  t 

and  setting  the  time  derivative  equal  to  zero  yields 

t  -  -(  +  <('3;S  -  t* 

The  corresponding  range,  the  "time-corrected  range,"  is 

1 

=  R(a,t*)  =  [(aj^t*  +  a2-x^*)^  +  (a3t*  +  84-7^*)^]^ 

The  asterisk  designates  values  corresponding  to  t  =  c*. 

Since  is  available  as  a  by  product  of  the  straight-forward 
CHURN  process,  the  additional  computation  needed  for  time-corrected 
range  is  trivial.  The  accuracy  of  this  range  is  considerably  better 
than  for  range  at  solution  time.  Since  the  time-corrected  range  is 
older,  however,  and  since  the  other  solution  parameters  are  not 
improved,  its  usefulness  is  controversial. 
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III.  COVARIANCE  OF  E 

As  previously  stated,  the  Bossard  method  makes  the  assumption 

* 

that  the  time  derivative  E  of  the  residual  E  is  exponentially  auto- 
correlated.  According  to  the  target  motion  used  by  Bossard  in 
reference  2,  zig  times  are  to  be  randomly  dravm  from  a  Poisson 
distribution  such  that  the  probability  of  no  zig  within  time  interval 


t£  to  tj  is 


-|ti-tj  l/t^ 


(3-1) 


where  tjjj  is  the  mean  time  between  zigs.  The  new  course  is  to  be 
randomly  selected  from  a  rectangular  distribution  extending  ±  6Cj  max 
degrees  about  a  mean  course.  With  these  conditions,  Appendix  B  shows 
that  the  target  course  deviation  has  a  covariance  matrix 

t 

2  j  “  -  tj/tn,  2 

€  (.  6  Cj  i  C,j,^)  £ j  =  e  (degrees)  (3-2) 


2  i  t  2  • 

where  "  3  ^  bearing  rate  is  small,  E  is 

nearly  proportional  to  4c.p.  Thus  the  covariance  matrix  for  E  should 

have  the  same  form; 


E  a  /•  ‘Tx  E  " 
C^  5  «(ee")  =  C“o  e 


(3-3) 


In  order  to  use  equivalent  weights  in  the  CHURN,  which  minimizes 

T  ’T  • 

E  WE  rather  than  E  WE,  the  corresponding  covariance  of  E  is  needed. 
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This  has  been  derived  as  outlined  here  (given  in  more  detail  in 
Appendix  C); 

Gy  ^  €(EiEj)  =  «  |^(Eo  +1  E  dt)(Eo  +  J  ^  E  dt)] 


(3-4) 


=  € 


[eq^  +  Eo(Ei  +  Ej  -  2Eo)]  +1^1^  <[E(t)E(t')]  dt  dt' 

2  rH  r*'j 


"^00  ^io  ^oj  *^g 


^  rr 


dt  dt' 


(3-5) 


After  integration  one  obtains  as  final  result 

+  /T  ^  ^  -  Iti  -  til/tm 

°ij  “  ’^00  +  +  C^-i  +  C 


+  e 


3.0  “OJ 


2  2  r 

E  L 

+  e  ^  -  1  +  (It^l  +  Itjl  -  |t^  -  tjD/tj^J 


(3-6) 


Tlie  functions  C^Q(t£),  CQj(tj),  and  the  constant  which  result 
from  integration  are  completely  arbitrary  as  far  as  the  original 

♦  C* 

specification  of  E  is  concerned,  but  must  be  selected  so  that  C 

retains  the  general  properties  of  covariance/autocorrelation  matrices: 

F  E 

C  should  be  positive  definite  as  well  as  symmetric.  Therefore 

and  Cqj  should  be  identical  functions  of  t^  and  t j ,  respectively. 

For  simplicity,  these  restraints  were  met  by  arbitrarily  setting 


_  _E  _  „E  _  rf  2  ^  2 

Cio  -  Cjo  -  Cqo  -  0^ 
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2  2 

Again  for  simplicity,  the  factor  0^  t^  was  taken  as  unity,  since 
for  weighting  purposes  a  constant  factor  on  the  covariance  has  no 
effect.  The  result  is  the  expression  used  'n  many  of  the  tests: 


C  Simple  *  -e  +  e 


m 


"  1 1  j  I  /  tjjj 

+  e  +  (Itj^l  +  Itjl  -  |t^-tj|)/tn, 


(3-7) 


The  origin  of  time  was  taken  at  the  middle  of  the  data  window. 

Seeking  a  more  logical  way  of  selecting  the  arbitrary  functions 

EE  E  ■~E 

Cio,  Gqj  and  constant  Cqq,  we  have  derived  a  covariance  of 

deviations  from  the  mean  E  over  any  selected  time  interval  t^  -  T  to 

tg  +  'T  .  The  derivation  (given  in  full  in  Appendix  D)  starts  with 


t^+T 

Ei 

a-^ 


dti  =  0 


(3-8) 


It  follows  that 


ft^+r 

<(EjE)  =  (1/2T)  J  (Ej^Ej)dtj^  =  0 


(3-9) 


0  =  (1/27)  /  C^j  dtj^  (3-10) 

•'t  -T 
‘■a  ' 

After  substituting  from  equation  (3-6)  into  (3-10)  and  integrating, 

E  E  E 

it  is  possible  to  solve  for  C^q,  and 
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+  (tj-t^)2/2t^T 


+  (tm/'O  +  'T/Stn,  +  /2t2)  1  -  e 


(3-11) 


Some  computer  runs  were  made  with  this  covariance  formula,  using  for  tg 
the  center  of  the  window,  and  for  T  one-half  the  window  length* 

There  is  some  logic  in  the  notion  of  placing  tg  at  the  end  of  the 
data  window  instead  of  at  the  center.  This  should  have  the  effect  of 
weighting  the  more  recent  observations  more  heavily,  possibly  yielding 
more  up  to  date  solutions.  This  idea  has  not  yet  been  tried. 

We  were  somewhat  puzzled  to  find  that  certain  widely-differing 
covariance  matrices  yielded  near  identical  range  results,  while  in 
other  cases  a  diminutive  change  to  the  assumed  covariance  matrix 
significantly  influenced  the  results.  The  former  situation  is  illus¬ 
trated  by  c|j  Simple  and j  for  which  equations  are  given  above, 

(3-7)  and  (3-11) .  Numerical  elements  of  the  covariance  matrices 
discussed  here  and  of  the  weighting  matrices  which  are  their  inverses 
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are  given  in  Appendix  E.  No  resemblance  is  apparent  either  between 
covariance  matrices  or  weighting  matrices,  yst  the  results  obtained 
were  nearly  identical.  The  other  situation  is  illustrated  by 
(Cy  Simple  +  .0021) .  Although  the  diagonal  elements  of  Simple 
are  thus  increased  by  no  more  than  0.2%,  the  corresponding  inverse  matrix 
is  thereby  noticeably  changed,  and  range  results  were  appreciably 
influenced . 

Experiment  has  also  shown  that  the  combination  (C^^j  +  .0021) 

E 

produced  ranges  nearly  identical  to  those  of  Simple  +  .0021)  , 

EMS  values  differing  less  than  5  yards  in  cases  tested. 

IV,  TEST  CASES 

The  tests  here  reported  used  mainly  the  geometry  shown  in 
Figure  2,  this  being  the  same  as  used  by  Bossard  in  reference  2. 

According  to  the  reference,  the  target  path  in  this  geometry  was 
constructed  by  the  algoritlim  described  in  section  III,  using  tjjj  =  15 
minutes  as  the  mean  time  betweei.  zigs,  and  selecting  each  new  course 
from  a  rectangular  distribution  lying  within  ±40°  of  the  mean  course, 
142°.  Own  ship  retained  an  unchanged  schedule  of  zigs  throughout  all 
tests,  and  both  ships  maintained  speed  constant  at  approximately 
10  knots. 

Solutions  were  computed  for  a  moving  data  span  or  "window" 
usually  20  minutes  in  length,  which  included  11  bearing  samples  taken 
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at  odd-numbered  minutes.  Thus  the  first  solution  came  at  21  minutes, 
followed  by  new  solutions  at  2-minute  intervals. 

In  order  to  test  certain  hypotheses,  variations  in  the  above 
conditions  and  solution  methods  were  introduced.  These  included  the 
straight  target  path  and  sinuous  target  path  shown  in  Figure  3,  bearings 
with  noise,  and  a  longer  data  window. 

V.  DISCUSSION  OF  RESULTS 
1.  Churn  vs .  Quadruples 

The  initial  question  attacked  was  why  the  quadruple  method 
proposed  in  reference  2  yielded  better  ranges  than  the  CHURN  for  the 
synthetic  case  used  for  demonstration  (zigging  target,  noise-free 
bearings) .  Several  ideas  were  tried  before  equal  performance  was 
achieved . 

In  the  quadruple  method,  the  single  variable  range  is  optimized. 

In  CHURN,  on  the  other  hand,  the  "time-corrected  range"  corresponds 
to  the  minimum  sum  of  variances  of  the  two  variables  x  and  y.  It 
seemed  possible  that  if  by  rotation  of  coordinates  the  range  vector 
was  made  nearly  parallel  to  the  y-axis,  and  if  then  the  variance  of  y 
alone  was  minimized,  that  a  better  range  value  would  be  obtained.  No 
significant  gain  was  realized,  however,  as  can  be  seen  by  comparing 
columns  3  and  4,  first  row,  of  Table  1. 
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The  next  experiment  was  intended  to  try  in  the  CHURN  a  weighting 
matrix  equivalent  to  that  used  in  the  quadruple  method.  Through  an 
initial  misunderstanding  of  this  weighting,  weights  corresponding  to 

Cij  =  e 

were  applied  to  the  CHURN.  Table  1  exhibits  the  poor  results  thus 
obtained  (column  5) . 

Another  experiment  used  the  same  weights  divided  by  the  range. 

It  is  well  known  that  CHURN  produces  biased  estimates  because  the 
residual  E  is  scaled  by  the  range.  In  a  synthetic  problem  in  which 
true  ranges  are  known  a  priori,  some  of  the  bias  can  be  removed  by 
dividing  out  the  range.  Column  6  of  Table  1  shows  the  improvement  to 
be  neglibible. 

Further  examination  of  Bossard's  quadruple  formulation  revealed 
that  the  assumed  exponential  (Poisson)  distribution  applies  essentially 

4 

to  E,  rather  than  to  E.,  The  derivations  of  section  III  and  Appendices  B, 
G  yielded  the  corresponding  covariance  matrix  needed  to  apply  this 
assumption  in  the  CHURN.  As  shown  in  the  first  row,  last  column  of 
Table  1,  the  range  errors  thus  obtained  with  the  CHURN  are  even  slightly 
smaller  than  those  of  the  quadruple  method.  It  was  concluded  that  the 
merit  observed  in  the  quadruple  method  arises  frcr.n  the  weighting,  rather 
than  from  geometrical  properties.  Subsequent  tests,  using  the  tHIURN, 
were  directed  toward  learning  the  effect  of  this  weighting  and  its 
variations  under  differing  conditions. 
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Although  achieving  a  substantial  advantage  for  a  maneuvering 
target  in  the  absence  of  bearing  meansurement  noise,  the  type  of 
weighting  used  in  the  quadruple  method  performed  poorly  when  noise 
v/as  present.  In  order  to  separate  the  effect  of  bearing  noise  from 
the  effect  of  maneuver,  rests  were  performed  on  a  straight-running 
target  with  unweighted  CHURN  and  with  weighted  CHURN.  The  second 
row  in  Table  1  shows  that  the  unweighted  CHURN  performs  better  for 
this  case. 
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2.  Var5.atioR  of  Covariance  Parameters 

We  desired  to  .  vestigate  whether  adjusting  the  parameters 

of  the  assumed  covariance  upon  which  x/eighting  is  based  would  bring 

about  any  improvement,  especially  in  the  presence  of  bearing  noise. 

The  second  line  Ln  Table  1  indicates  that  the  Bossard  weighting 

performs  poorly  xAen  uncurrelated  bearing  noise  is  present.  It 

seemed  reasonable  that  other  values  of  t^j  might  work  at  least  as 

well,  since  only  three  zigs  at  most  occurred  within  any  window,  and 

such  a  zig  pattern  could  arise  from  a  wide  range  of  tm  values  with 

almost  equal  probabilities.  In  addition,  the  derivation  process 

(Appendix  C)  contributes  to  the  covariance  of  E  an  undetermined 

constant  Cqq  and  undetermined  functions  C^^q,  Cj^.  A  series  of  tests 

was  performed  to  compare  results  with  variations  in  these  parameters, 

and  the  range  errors  are  given  in  Table  2.  The  upper  half  of  the  table 

E  £  E 

shows  results  obtained  by  changing  (C^^  +  C^j  -  C^q) ,  e.g.  adding  a 
uniform  constant  to  each  element  in  the  covariance  matrix.  The 
effect  on  range  errors  is  insignificant,  with  or  without  noise. 

The  lower  part  of  Table  2  indicates  the  effect  of  varying  tj^^. 
Performance  in  the  no-noise  case  deteriorates,  but  in  the  more  realistic 
case  of  bearings  x?ith  noise,  very  substantial  improvement  can  be 
obtained  by  reducing  tj^  to  0.2  -  0.1.  The  logic  underlying  this 
result  is  only  partly  clear.  The  shorter  mean  time  to  zig  implies  less 
correlation  between  observations,  and  thus  corresponds  better  to  the 
uncorrelated  bearing  noise  portion  of  the  residual  E. 
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TABLE  2 

VARIATION  OF  Al^D  (0?^,  + 


E  E  E 

Cio  Coj  -  Coo 

t^jCmin) 

No 

Noise 

Noise 

II 

o 

» 

O 

At 

Solution 

Time 

Time- 

Corrected 

At 

Solution 

Time 

Time- 

Corrected 

0.2 

15.0 

1914 

699 

6021 

6715 

1.0 

tf 

1915 

700 

6013 

6715 

5.0 

It 

1921 

709 

6040 

6716 

i.O 

28.6 

2003 

795 

II 

15.0 

1915 

700 

6013 

6715 

It 

7.0 

1995 

725 

II 

2.18 

4140 

3471 

n 

1.66 

3786 

2780 

u 

1.24 

3614 

2288 

II 

0.87 

4595 

2564 

3691 

2104 

II 

0.87 

4263 

1963  (10) 

II 

0.67 

3849 

2084 

M 

0.43 

4612 

2348  (10) 

Target  maneuver:  zigs 
C  =  -e 


+  dtj  +  Itjl  -  l)/tn 


^  ^10  ^  ^Oj  ^00 

(10)  RMS  of  10  noise  samples.  Other  entries  RMS  of  3  noise  samples. 
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3.  Covariance  for  Combined  Error  Sources 

Derivation  of  the  covariance  as  for  the  results  of  Table  2 
accounts  only  for  residuals  arising  from  zigs,  whereas  the  covariance 
commonly  used  in  the  CHURN  (proportional  to  an  identity  matrix) 
assumes  residuals  arising  only  from  Gaussian  bearing  measurement 
noise.  When  both  sources  contribiite  to  the  residuals,  it  is  reasonable 
to  add  their  covariance  contributions.  The  sum,  expressed  in  a  form 
to  make  the  units  of  both  contributions  consistent,  is 

2  2  E  0 

Combined  covariance  =  O'*  t  C  Simple  +  O  ~  i 

Em  EB 

where  C^  Simple  is  defined  by  equation  (3.7) ,  ffgg  is  the  variance  of 
E  caused  by  bearing  measurement  noise,  and  I  is  the  identity  matrix 
of  the  same  order  as  C^.  Rearranging, 

2  9  T?  9  2  2 

(Combined  covariance) /(  cr*  t  )  =  C  Simple  +  )I 

b  in  Eg  b  in 

from  which  wo  define  the  equivalent  terms 

C^  Combined  =  C^  Simple  +  Ol 

We  have  performed  a  number  of  CHURN  tests  to  determine  whether  the 
C^Combined  leads  to  significant  improvement.  Since  the  factors  upon 
which  ct  depends  would  not  be  known  with  precision  in  practical 
situations,  we  also  wished  to  determine  whether  its  value  is  critical. 
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Figures  4,  5,  6  show  the  effects  of  variation  in  ct  on  range 
errors  for  three  types  of  target  motion.  All  cases  include  uncorrelated 
bearing  measurement  noise.  The  RMS  range  errors  have  been  plotted 
against  in  order  to  normalize  for  the  several  values  of  tjj^ 

included.  Figure  4  verifies  an  expected  result,  that  for  a  straight¬ 
running  target  the  unmodified  identity  weighting  (  «  =  oo)  is 
optimum. 

Figure  5  applies  to  the  zigging  target  used  in  reference  (2) . 

This  graph  confirms  that  a  proper  blend  of  Simple  and  I  is 
significantly  better  than  either  alone  for  this  target  in  the  presence 
of  bearing  noise,  except  when  tj^  is  very  small.  It  also  indicates 
that  ct  may  vary  by  a  factor  of  three  either  way  from  the  optimum 
without  serious  degradation  of  the  range  solution. 

Figure  6,  for  a  sinuous  target  motion  and  bearing  noise,  also 
shows  advantage  for  the  Combined,  even  though  the  statistical 
distribution  of  the  residual  E  generated  by  this  maneuver  is  entirely 
different  than  that  generated  by  a  zigging  target.  The  factor  c*  is 
again  noncritical. 

Returning  to  consider  Figure  4  again  along  with  Table  2,  it 

p 

appears  that  as  much  is  gained  by  reducing  t^^^  as  by  using  C  Combined. 
Figure  7,  however,  supports  the  use  in  the  absence  of  noise  of  a  value 
near  the  theoretical  one,  tj^  =  15  minutes.  The  conclusion  to  be  drawn 
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is ,  that  the  theoretically  appropriate  values  of  tj^  and  a  work  at 
least  as  well  as  any,  but  when  the  effect  of  random  noise  is  large, 
certain  other  values  also  work  as  well.  A  more  specific  characteriza¬ 
tion  at  this  stage  seems  hazardous. 

4.  Tactical  Selection  of  Covariance  Parameters 

It  is  clear  from  Figures  4  through  7  that  for  best  results 
from  the  CHURN,  the  assumed  covariance  parameters  tj^  and  d  should 
be  adjusted  to  suit  the  noise  and  target  maneuver.  Since  the  a  priori 
information  on  target  maneuver  available  to  the  tracking  ship  would 
be  limited,  a  method  of  adjustment  based  on  preliminary  data  analysis 
is  desirable.  We  have  examined  several  schemes  for  selecting  the 
parameters,  with  doubtful  or  negative  results.  One  method  consisted 
of  comparing  the  variance  of  bearing  residuals  based  on  the  same  data 
but  with  different  assumed  functions  for  the  covariance  of  E.  It  was 
hoped  that  small  residuals  would  correlate  closely  with  small  range 
errors.  In  Figure  8,  the  variance  of  residuals  is  plotted  against 
RMS  error  of  time-corrected  range,  and  in  Figure  9  against  RMS  error 
of  soluticn-time  range.  Points  plotted  in  these  figures  are  coded  to 
distinguish  three  combinations  of  zig  and/or  bearing  noise,  and  for 
each  combination  a  variety  of  covariance  functions  is  represented. 
While  some  correlation  is  noticeable  between  bearing  residuals  and 
range  errors,  the  relation  is  not  as  consistent  as  one  desires  as 
a  basis  for*  choice. 
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5.  Zig  Detection  by  Autocorrelation  of  Bearing  Residuals 

Unfavorable  results  were  obtained  from  our  attempts  to  detect 
zigs  by  means  of  the  autocorrelation  of  bearing  residuals.  The 
principle  involved  is  that  the  shape  of  the  autocorrelation  curve 
estimated  from  the  sample, 

R(T)3jjjj,pig  =  [mean  jB(t)  6B(t  +  T)]/[mean  dB(t)  5B(t)] 
should  be  different  for  zigging  and  nonzigging  targets.  If  one  assumes 
that  bearing  residuals  should  have  the  same  covariance  (except  for 
a  scale  factor)  as  the  cross-range  residual  E,  then  for  a  zigging 
target  with  bearing  noise  the  expected  R( I") sample  would  be  the 
theoretical  curve  drawn  in  Figure  10  on  the  left  side.  This  curve 
was  computed  from  the  covariance 

+  .002  I 

For  a  non-maneuvering  target  with  noise,  the  expected  value  of 
R(t)  sample  would  be  unity  at  the  origin  and  zero  every  where  else. 

Figure  10  exhibits  the  empirical  R(T)ga^p]^g  points  obtained 
for  two  sets  of  three  runs.  Both  sets  are  identical  in  conditions 
except  that  one  has  a  zigging  target  while  the  other  has  a  straight¬ 
running  target.  The  "random"  sequence  of  bearing  noise  values  is  the 
same  for  each  set.  It  does  not  appear  that  any  clearcut  characteristic 
distinguishes  the  zig  cases  from  the  straight  path  cases.  Furthermore, 
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each  curve  is  an  average  over  49  minutes  of  data,  which  is  generally 
too  long  to  wait  for  the  basis  of  a  shipboard  decision.  \  ith  a  lower 
noise  level  a  better  result  would  possibly  be  obtained,  but  the  level 
assumed  is  not  unreasonable.  The  value  =  0.1°  applies  to  generated 
bearings  two  minutes  apart,  corresponding  to  a  mean  of  sixty  or  so  raw 
bearings . 

6.  Time -Corrected  Solutions  vs.  Longer  Windows 

Reference  (2)  proposes  extrapolation  from  a  sequence  of  time- 
corrected  ranges  to  obtain  range  at  present  (or  desired)  time.  It 
gives  examples  of  time-correction  applied  to  CHURN  to  support  the 
assertion  that  the  line  joining  two  time-corrected  ranges  from  suitably 
chosen  windows  can  be  extended  to  estimate  present  range  at  least  as 
well  as  a  single  solution  covering  the  entire  interval.  The  advantage 
is  said  to  be  greater  if  the  tracking  ship  executes  a  lag-lead-lag 
series  of  course  legs. 

Our  tests  of  this  method  employed  different  conditions  than 
those  ul  reference  (2) .  The  most  significant  differences  are  that  we 
used  the  zigging  target  of  Figure  2  (instead  of  a  straight  path)  ,  a 
different  noise  level,  and  weights  corresponding  to 
C^Combined  =  C^Simple  +  .002  I 

in  the  CHURN  (instead  of  the  unweighted  CHURN).  This  covariance  has, 
in  our  other  tests,  yielded  the  best  overall  range  results  for  this 
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geometry  and  noise  level  (  Cg  =  0.1°),  Time-corrected  ranges  were 
obtained  for  every  20-minute  window  possible  in  the  run,  and  solution 
time  ranges  for  every  28-minute  run.  The  process  was  repeated  for  two 
additional  runs  with  different  noise  samples.  In  Figure  11  the  solid 
lines  represent  errors  of  time-corrected  range  plotted  against  the 
"best  range  time"  at  which  the  time-corrected  range  applies;  certain 
points  have  been  labeled  with  the  time  of  solution. 

As  predicted  by  reference  (2) ,  the  best  range  time  tends  to  fall 
near  the  beginning  of  a  window  which  spans  a  lag-lead  sequence  of 
tracking  ship  legs  (points  labeled  21  or  39) ,  and  near  the  end  of  a 
window  which  spans  a  lead-lag  sequence  (29  or  49) .  Thus  the  pairs 
21  and  29  or  39  and  49  provide  favorably  long  baselines  for  extrapola¬ 
tion.  Circled  points  terminating  the  short-dashed  lines  show  the 
range  errors  resulting  from  extrapolation  to  present  time,  at  29  or  49 
minutes , 

Figure  11  also  shows  the  range  errors  at  solution  time  using 
28-minute  windows  (square  points) .  The  28-minute  windows  ending  at 
29  or  49  minutes  encompass  nearly  the  same  data  as  the  pairs  of 
20-minute  windows  used  in  extrapolation.  At  29  and  49  minutes, 
therefore,  a  direct  comparison  can  be  made  between  the  two  methods: 
extrapolation  of  time-corrected  range  using  pairs  of  short  windows, 
versus  solution  time  range  using  longer  windows.  In  this  case  the 
verdict  is  nearly  a  tie,  but  it  is  evident  that  the  scatter  caused 
by  bearing  noise  would  mask  the  small  advantage  which  might  exist 
for  either  method. 
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VI.  CONCLUSIONS 

We  find  that  the  major  contribution  of  Dr.  Bossard  in  reference  (2) 
to  passive  TMA  technology  is  to  demonstrate  the  effectiveness  of 
appropriate  statistical  weighting.  Our  own  tests  have  shown  that  the 
CHURN,  with  equivalent  weighting,  performs  as  well  as  the  quadruple 
method  presorted  in  reference  (2).  The  importance  of  weighting  would 
apply  to  any  type  of  solution  using  redundant  data. 

The  weighting  used  in  reference  (2)  was  derived  to  fit  a  particular 
target  maneuver  strategy,  assuming  no  bearing  noise.  When  random 
bearing  noise  is  present,  however,  solutions  with  this  weighting  are 
poor  even  without  zigs;  indeed,  the  unweighted  version  of  CHURN  is 
then  optimum.  A  combined  weighting  has  been  found  to  be  best  when 
both  noise  and  maneuver  are  present.  Although  the  combining  proportions 
are  not  critical,  further  work  is  needed  to  develop  a  method  for 
selecting  the  best  weighting  on  the  basis  of  information  available 
to  the  tracking  ship. 

Compared  to  the  unweighted  CHURN,  the  CHUkN  with  weighting  derived 
for  a  zigging  ta^^et  has  been  found  superior  for  a  sinuous  target  path, 
although  the  residuals  generated  by  sinuous  motion  have  an  entirely 
different  statistical  distribution.  ^  This  result  leads  us  to  hope 
that  one  type  of  weighting  can  be  used  for  a  variety  of  target 
strategies. 

Our  attempts  to  distinguish  target  zigging  from  straight  path 
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by  examining  bearing  error  autocorrelation  were  unsuccessful-  This 
idea  would  probably  worK  at  lower  levels  of  random  bearing  noise. 

In  our  tests  of  time -corrected  ranges  extrapolated  to  estimate 
present  range,  no  advantage  was  observed  over  solution  time  range. 

Since  the  number  of  cases  tried  was  only  three,  the  Scatter  due  to 
noise  masked  whatever  small  difference  might  exist.  Time-corrected 
ranges  probsbly  have  other  applications,  since  they  indeed  have 
smaller  variance  than  solution-time  ranges.  For  instance,  information 
that  the  target  has  been  approaching  since  best  range  time  may  exist 
in  a  form  not  available  to  the  computer;  even  if  solution  time  range 
is  wild,  the  time-corrected  range  would  provide  a  fairly  reliable 
upper  limit. 

The  statistics  of  maneuver  introduced  here  in  the  CHURN  can  also 
be  introduced  in  the  optimal  (Kalman)  filter  technique  now  coming  into 
favor.  Reference  3  describes  a  Kalman  filter  and  trials  on  real  data. 
The  formulation  given  permits  statistical  variations  of  target  velocity 
components,  i.e.  target  maneuver,  and  in  some  of  the  trials  the  filter 
did  assume  non-zero  statistics  for  velocity  changes.  All  of  the  cases 
tried,  however,  had  a  non -maneuvering  target.  The  suboptimal  results 
correspond  to  our  own  results  with  the  CHURN  (reported  herein)  ,  wlien 
using  weighting  appropriate  to  a  zigging  target  on  a  straight-running 
target. 
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VII.  AREAS  FOR  FUTURE  WORK 

Taking  full  advantage  of  weighting  requires  a  method  for 
selecting  a  near-optimum  weighting  function,  given  such  information 
as  a  tracking  ship  can  possess  in  the  tactical  environment.  Our  few 
unsuccessful  attempts  were  aimed  at  extracting  a  weight  selection 
criterion  from  bearing  data  alone.  Vfhile  we  believe  this  not  to  be 
hopeless,  further  studies  should  consider  the  use  of  bearings  together 
with  other  sensor  data  in  a  combined  solution  method.  With  the 
maneuver  detection  problem  thus  greatly  reduced,  appropriate  bearing 
statistics  would  be  available  to  enhance  the  combined  solution. 

In  addition  to  target  zigs,  other  sources  of  correlated  errors 
exist  for  which  appropriate  weighting  could  improve  solutions.  Examples 
are  own  ship  position  keeping,  and  bearing  measurement  bias.  Enough 
is  known  about  these  particular  sources  so  that  error  modeling  would 
not  be  difficult.  Investigation  is  needed  to  determine  how  large  the 
potential  improvement  would  be,  and  how  to  introduce  the  statistics 
into  a  solution. 

Inasmuch  as  the  Kalman  filter  solution  method  appears  soon  to 
become  the  standard,  further  innovations  should  either  be  implemented 
in  this  framework  or  tested  in  competition  to  it  with  respect  to 
accuracy  and  computing  economy. 
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APPENDIX  A 

Quadruple-Bearing  Range  Solution 

The  following  is  a  computing  algorithm  for  obtaining  range 
estimates  from  sets  of  four  observed  bearings,  and  for  averaging  with 
appropriate  weights  the  estimates  spanning  a  data  window.  The  result 
is  a  range  solution  at  best  range  time,  i.e.  a  time-corrected  range. 
While  different  in  form,  the  equations  used  here  are  equivalent  to  those 
presented  by  D.  C.  Bossard  in  reference  2. 

In  a  time  interval  (window)  containing  LW  discrete  equally-spaced 
observations,  the  observation  times  are  designated 
tjj  =  time  at  beginning  of  window 

H  =  tb+i  +  i4t 

Similar  subscripts  on  other  parameters  refer  to  corresponding  observa¬ 
tion  times.  Further  definitions  are 
B^  =  bearing 

Xoi,  =  own  ship  coordinates 
S^j  =  sin  B^  cos  Bj  ••  cos  B^  sin  Bj 
G..  =  cos  Bj  cos  Bi  +  sin  B.  sin  B^ 

Dij  =  (xo.i-xoj)cos  B.  -  (y^.-yoj)sin  Bi 
hj  =  -C^oi'^oj^®^"  h  -  (yoi-yoj)<^os  Bi 
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Using  the  above  definitions  compute  the  following  for  1  i  (LW-2): 


®i  ~  ^i,i-l  "  ^i-l,b  ^i,bfl  ^W-l,b 

h  •  '■h-l  -  l^b) 

’'l  ■  “l-l.b  Si,b+l/% 


The  range  estimate  for  a  single  set  of  four  bearings,  although 
not  explicitly  used  in  the  succeeding  computation,  is  cos  £  j^) 
in  the  equation 

(^bCOs  <b)i=  «i+  ^iCST^b+l 

^■^i-l,b‘^i,b+l^b-b3,b’  ^i-l,b‘^i,b-l»  ^i,b»  ■^i-l,b^ 


‘"b+1 

^i-l 


•r  Di 


in  which  S^I  is  the  line  of  sight  component  of  target  speed,  and 
is  the  angle  subtended  by  the  cross  range  residual  Ejj. 

We  proceed  next  with  the  computation  of  weight.  It  is  shown  in 
reference  (2)  that  the  in-line  range  residual  is  approximately 
proportional  to 

jRbi  (-Eb  +  Ei,.p/D^ 
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This,  together  with  the  assumption  that  the  covariance  of  E  is 


-Iti'tj  l/t^  ^ 
e  =  K 


leads  to  a  covariance  for  the  in-line  range  residual 


R  ii+ll  Ij+ll  li-jl 

Ci  j  =  (1-P  -P  +P  )/(D^Dj) 


This  is  a  symmetrical  matrix  in  idiich 

1  i,j  sj:  (LW-2) 

Compute  the  matrix  and  its  inverse.  Compute  weights  using  the  formula 


LW-2 

-1 

/LW-2 

LW-2  -1 

r 

(C^)  1 

/  L 

Z  (c^) 

ij/ 

i=l 

ij 

Compute  the  weighted  means 


_  LW-2 
Cl  =  £ 

i=l 


_  LW-2 

?  =  ,£  "i  Pi 
i=l 

_  I^W-2 

7’  =  r  w.  y. 
i=l 


t  =  tj,  +  p 
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Interpolate  for  B  =  3(t) .  Then  compute 

t*  =  t  +  y  cos(B-Bjy{.TL) /cos  (B-Bb) 
luCerpoiate  for  3*  =  B(c*) ,  and  finally  compute 

^  ■  ^b+l,t* 

+  (  a  -  Ljj— )  cos(B  -  ®b+l^  /cos(B-Bb)]  /cos(B*-Bb4.i^) 
which  is  the  time-corrected  range  for  the  window. 
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APPENDIX  B 

Justification  of  Exponential  Correlation  for  E 


Consider  a  target  vessel  which  is  traveling  at  constant  speed, 
roalcing  good  a  mean  course  while  performing  zigs  as  evasive  maneuvers. 
The  zig  strategy  is  assumed  to  be  as  follows:  the  time  for  each  zig 
is  selected  so  that  the  probability  of  no  zig  in  die  time  interval 


tj^  to  tj  is 


- l/tm 


P(no  zig  between  t^  and  tj)  =  e  (1) 

and  the  deviation  from  mean  course  resulting  from  each  zig  is  drawn 

2  1 

randomly  from  a  rectangular  distribution  of  variance  ‘^CT  3  (  *  max) 
Then  the  covariance  of  target  course  deviation  6c^  is 


C^j  =  ^  [6c.j>(t£)  iC.j(tj)] 

2 

“  between  t^  and  tj) 

+  (zoro)P(any  zig  between  t^  and  tj) 
2  -iH-tjl/tm 

=  Oct  e 


(2) 

(3) 

(^) 


The  constant  t^  turns  out  to  be  the  mean  time  between  zigs. 

As  a  linearized  approximation,  the  cross -range  error  residual  E 
resulting  from  such  target  maneuvers  may  be  written 


^i  “  [  ^TI 


(5) 
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APPENDIX  C 

* 

Covariance  of  E  Derived  from  Covariance  of  E 

Given  that  the  covariance  of  the  time  derivative  E  of  the  error 
residual  E  is 


e..  =  o~  e 
13  E 


l/t„ 

i  J  Ul 


(1) 


then  the  covariance  of  the  residual  itself  is 


Cfj  =  €(EiEj)  =  tj^CEo  +J  Edt)(E^,  +1  ^  Edt)j  (2) 

=  <[eo^  +  Eo(Ei  +  Ej  -  2Eo)J  ^  €[^E(c)E(t ')  ]  dtdt'  (3) 


dt  dt' 


Integration  will  be  performed  piecewise  over  regions  within  which 
the  sign  of  the  exponent  can  be  insured.  The  double  integral  becomes, 
if  Cj  0 

|tjj-|t’  e(t-t')/tn,  ^t  ^|‘■i 

tj  r  -f/t^ 


dt' 


2  ,  -4/tm  -tj/tm  , 

=  tm  (-e  +  e  +  e  -  1)  +  2trotj 


(5) 
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If  tj  ^  ^  0,  symmetry  with  respect  to  i  and  j  leads  to  the 

analog  of  (5) 

-(tj-ti)/tjjj  -tj/tjg 

tj2  (-e  +  e  -  1)  +  2ti„ti 


(6) 


If  0  Cj  y.  t^,  the  integral  of  (4)  becomes 


ro  r  ^t'  (t-t')/t^  ^o  -(t-t')/tjj, 

e  dt  +  j  e 

4' 


^t.  L-'ti 


dt 


dt' 


Of  (ti-t')/t^  f/t^ 


0- 


-  e  +2 


dt' 


tj^  (-e  +  e  +  e  -  1)  -  2tmtj 


(7) 


By  symmetry,  if  0  ^  tj:  "?/  tj  the  analog  of  (7)  is 

-  (tj-ti)/tjj,  tj[/t[jj  tj/tjj, 

t^  (-e  +  e  +  c  -  1)  -  2t^tj^ 


(8) 


If  tj^  ^  0  ^  t j ,  the  integral  of  (4)  becomes 


rtjT  rt-i  -(t-c')/tjn 

■o  *-4 


dt 


dt' 


.tj  f  -(ti-t')/t„  t'/tjn-i 

H  r 


dt' 


t^  (-e  +  e  +  e  -  1) 


(9) 
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APPENDIX  D 

Covariance  of  E  for  E  =  0  Over  Selected  Period 


In  a  previous  derivation  the  expression 


+  cr. 


V 


-  It.-t.  l/t 
-e  ^  ^  + 


-itj/t, 


m 


+  e 


+  (Itil  +  Itjl  -  ft) 

was  obtained  for  the  covariance  of  E,  corresponding  to  a  given 
covariance  for  the  tine  derivative  of  E: 


2  -i'=r'=ji/Su 

<7  •  e 


(2) 


E  E 

C^Q  stands  for  a  function  of  tj  only,  Cqj  a  function  of  tj  only,  and 

E 

Cqq  a  constant.  These  arbitrary  quantities  are  not  further  restricted 

by  equation  (2) ,  nor  by  the  general  requirements  for  autocovariance 

£ 

functions  e.xcept  that  C^^j  must  bo  positive  definite.  In  order  to  get 
explicit  expressions  for  the  arbitrary  quantities,  another  requirement 
which  seems  reasonable  is  here  imposed:  that  the  mean  value  of  E 
over  a  selected  interval  t_  -  T  to  t  +  T  is  zero, 

o  a  » 


t +r 


■-  f 

|J-T 


Ej^dtj  ”  0 
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The  covariance  which  meets  this  requirement  will  be  designated 'c®. 
From  (3)  it  follows  that 

1 

€(EjE)=~y  £(E.Ej)dt.  =0  (4) 


t„-i-T 


rv  1  r  a'  ■ 

^  =  C.  .dc. 

2t  i  LJ  1 


■'t  -r 
a 


Before  substituting  (1)  into  (5) ,  equation  (1)  is  abbreviated  by 
replacing  with  F^,  Fj ,  and  K  all  of  the  terms  in  t^  only,  terms  in 
Uj  only,  and  constant  terms,  respectively.  Equation  (1)  becomes 


2  2  [  ~  t  I  / 

E  V  [-^  +  +  +>^ 


which  is  substituted  into  5,  dropping  constant  factors; 

[-C  +  +  +Kjdti 


*^3  r 


+  (tj^-t 


rVrr  -(q-tp/t^ 

+  J  -a  -  (ti-tj)/tm 


>•  +  T 

.  .._T  I 


("'a 

+  I  F^dtj  +  2TFj  +  2TK 


^■ta-r 
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y 


■V  +  q/2t^  - 


■’‘j 


r  *1-1=1) /tm  , 

^  *'in®  “ 


nt^-t-r 


*^3 


(9) 


I  F.dt.  -t-  2tF.  +  irYi 
-'t  -T  J 


1-e 


(ta-T'-tj)/t„' 


2t_ 


m 


tj(t,-T) 


m 


-(t,+  T-tj)/t:„n 

1  9  0 

/  r  oi 

*Tn 

e 

-[(t,+  T) 

/2tm  +  [t.(t^+T)  -  tj 

F.dt.  +  2rF.  +  2rK 

*  *  J 


ta-'T 


(10) 


0  = 


tmi  e 


!  +  e  ' 


-  2 


+  ('=j  - 


(11) 


I  Fidti  +  2rF.  -f  2rK 
J 
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By  symmetry,  (11)  may  also  be  written  with  i  and  j  interchanged: 


0 


-  (ti-ta)^/t^  -  r^/t^ 


rta+T 

+  /  F.dt.  +  2TFi  +  2rK 
ta-T 


Subtracting  (11)  from  (12),  noting  that  jF^dtj^  =  /Fjdtj: 


-(ti-ta)^/ttn+  <t.-t^)2/t„+  2T[Fi-Fjj 


Because  t^  and  tj  are  independent,  the  terms  in  tj^  must  sum  to  zero  or 
a  constant,  and  those  in  tj  also  must  sum  to  zero  or  to  an  offsetting 
constant.  We  assume  that  the  constant  part,  if  any,  is  assigned  to 


K.  Then 

-'^/tmf  (ti-ta)/t„  -(ti"ta)/t„ 
e  e  +  e 


-(ti-ta)^/tj„  +  2TFj^  =  0  (14) 


,  M  ^  ^V‘=a)/‘4n  -(ti-ta)/tm1  , 

-C*m/2T)e  e  ^  +  e  +  (tj^-t^) /2t^T  (15) 


An  analogous  expression  applies  for  Fj. 
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To  evaluate  K,  substitute  (15)  and  its  j-analog  into  (11): 


0  =  W 


t  +  e 


m 


•}  -  <'j-' 


-2\-  (t,-tg,)2/t^  -T^/t^ 


r  ^  ■  I .  /o  ^ 

^*'i“*'a^^*'tn  “(*'i"*'a^ /*Tii 

e  +  e 

Vt  ^ 

+  (ti-tj2/2t^^|dt.  +  2T|-(t^/2T)e 


+  (tj-t^)V2t^T  }.  +  2TK 


-^/tm 


+  e 


(16) 


-e 


(ti-ta)/t^  -(ti-ta)/t„] 
+  e 


ta+T 


V'’* 


+  (1/t^T) 


tJ/6  -  tgt^/2  +  t^t^/2 


t^+T 


Jt^-T 


+  2TK 


(17) 


2  9  -'T/t 

0  =  -2t^  -  Tlt^  -  (4/t)  e  l^e  -  e 


m 


+  (l/6t^T)[(t3+T)^  -  (ta-T)^J  -  (ta/2t^r)J(t^+T)^  -  (tg-T)^! 

+  (ta/2t„r)[ 


(t^+T)  -  (t^-T)  j  +  2rK 


(18) 
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=  t^/r  +  r/3t^  +  (4/2  r^)  [l- 


-2T/t 


If  equation  (15)  and  its  j -analog  and  (19)  are  substituted  into  (6) , 
the  complete  expression  for  covariance  Is  obtained; 

— p  2  2  f 


-  (t^/2T)e  e 


e  +e  +  (ti-ta)^'2t^r 


-  it^/2r)e  e 


”^‘^“rar  “^'a^  ^*Tn  2; 

e  ^  +  e  ■’  J  +  it.-t^)^l2t^r 


2  2  r 

+  t^/T  +  T/3t^+  (4/2T  )[^1  -  e 


If  is  taken  as  the  time  origin  and  the  window  of  length  w 
coincides  with  the  interval  tg-T  to  tg+T  ,  then  (20)  simplifies  to 
V  0  n  f  -Iti-t. l/t 


(fin/w)  e 


-w/2t„,r  tj^  -t^  tj  -tj-j 

Q  +  e  +  e  +  e  J  +  (t^Z  +  t/)/t^w 

■  ] } 


+  2t^hi  I-  w/6t„j  + 


(2VV)  _ 


1  -  e 


Equation  (21)  has  been  used  for  some  of  the  reported  tests. 
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a 


1' 


^3 


a 


A 


a  (matrix  notation) 


ArFENDIX  ? 

Glossary  of  Symbols 

east  and  north  components,  respectively,  of 
target  ceLscity. 

east  and  north  comocnonts  resoectively,  of 
target  petition  with  re" 'act  to  an  arbitrary 
f2_:ed  oii-gin. 


^i' 

^2 

as 


A  (matrix  notation) 


B,- 


C  (Appendix  A) 


G^,  C^,  etc. 
Simple 


an  n-row  matrix  having  rows  of  the  form 
tj^  cos  cos  Bi,  -t^  sin  Bj^,  -sin  B^, 
observed  bearing  to  target  at  time  t^^. 
an  abbreviation  for  several  terms ,  defined  as 
introduced , 

covari.ance  of  the  superscribed  variable, 
a  particular  formula  for  the  covariance  of  E, 
derived  in  Append' 7:  C'e 

a  particular  fornn.la  for  the  covariance  of  E, 
derived  iv<  Appxnd'/:  D. 
target  couisi. 
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D  (Appendix  A) 


d  (matrix  notation) 


E 


g  (matrix  notation) 

I  (matrix  notation) 

L  (Appendix  A) 

P  (Appendix  A) 

R 

RVf 

R(T) 

'^)sample 

^TI 

S  (Appendix  A) 


an  abbreviation  for  several  terms ,  defined 
as  introduced. 

an  n-element  column  vector  with  elements  of 

form  x„,-  cos  B-  -  sin 

oi  1  ■'oi  oi 

the  cross  range  miss  distance  of  the  estimated 
bearing  line.  In  matrix  notation,  a  column 
vector  with  elements  E^. 

the  constant  term  in  the  matrix  representation 

of  CHURN  normal  equations. 

the  identity  matrix,  having  all  ones  on  the 

principal  diagonal,  all  zeros  elsewhere, 

an  abbreviation  for  several  terms,  defined  as 

introduced. 

an  abbreviation  for  several  terms,  defined  as 

introduced. 

range  to  target. 

estimated  range  at  best  range  time,  e.g,  time- 
corrected  range, 
autcacrrelation  function, 

estimate  of  R(T)  based  on  a  particular  sample, 
component  of  target  speed  along  bearing  line, 
an  abbreviation  for  several  terms,  defined 
as  introduced. 
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t 

t* 


W  (matrix  notation) 

Xq  .  Fo 

Ot 

a ,  0  ,  y  (Appendix  A)  - 

3 

€ 

€  jj  (Appendix  A) 

9  (Appendix  C) 

P 

2  2 

<T^  ,  ,  etc. 

T 

Ip  (matrix  notation)  - 


time. 

mean  time  interval  between  zigs. 
best  range  time:  time  at  which  the  expected 
range  error  of  a  solution  is  least, 
a  weighting  matrix,  usually  the.  reciprocal 
of  the  assumed  covariance  of  residuals, 
coordinates  of  tracking  ship  with  respect 
to  an  arbitrary  origin. 

a  combining  factor  used  when  adding  covariance 
contributions  from  target  zig  and  from  bearing 
noise. 

an  abbreviation  for  several  terms ,  defined  as 
introduced . 

(used  as  prefix)  indicates  error  quantity, 
expected  value  operator. 

angle  subtended  by  the  cross  range  residual  Ejj. 
symbol  for  an  unknown  expression,  dropped  when 
the  expression  is  determined. 

variance  of  the  parameter  indicated  by  subscript, 
time  interval. 

4x4  matrix,  coefficient  of  the  vector  a  in 
the  CHURN  normal  equations. 
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13.  A8STRAC*f 

12.  SPONSORING 

^.IILITARY  ACTIVITY 

(C)  In  a  new  sonar  bearings-only  solution  method,  Dr..i-.D..-G-v~BGSs*)f>d-4f^)e^^ 
'Wago«r--As'S064d4€S>  achieved  quite_sftac±'aGalar  reduction  in  range  errws  on  a  zigging 
target,  one-sixth  those  of  the  usual  (unweighted)  CHURN  method.  HHs  method  yields 
time-corrected  range  (value  at  time  when  expe^^  error  is  least)  and  weights  observa¬ 
tions  according  to  assumed  zig  statistics,  -i^^rd  also  advocates  extrapolating 
favorably-chosen  time-cprrected  ranges  to  obtain  present  range.  ,  -n., 

^  ^ 

We^i^  that  the  CHURN,  with  weights  equivalent  to  BosSjHls,  achieves  equally 
small  time-corrected  range  errors,  and  errors  at  solution  time  dne-third  those  of  the 
usual  CHURN.  Random  bearing  noise,  however,  seriously  degrades  solutions  using  Bos^hli' 
weights,  even  without  zigs,  in  which  case  the  unweighted  CHURN  is  optimum.  For  combi¬ 
nations  of  zigs  and  bearing  noise,  optimum  combined  weighting  functions  exist. 

y 

(p')  Unsuccessful  attempts  were  made  to  use  data  available  to  the  tracking  ship, 

|.g.  autocorrelation  of  solution  residuals,  for  selecting  optimum  weighting.  Auto- 
lorrelation  was  also  probed  for  zig  detection  clues  without  success. 

(.C)  Results  obtained  by  extrapolating  pairs  of  time-corrected  range  to  present 
time  were  about  equally  as  good  as  from  single  solutions  using  the  same  data. 

conclude,' that  B^oftrtH-s  important  contribution  “‘'s  to  show  the  effectiveness 
of  appropriate  statistical  weighting.. 
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